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\ A phenomenological estimate is derived such that the binding en- 

■ ergies of dimesons are expressed as combinations of masses of different 

r-j ■ mesons and baryons. The estimate is almost model-independent, the 

Qh! only major assumptions being that the wave functions of the two light 

O j' quarks in A c , A& and in the ccqq and bbqq dimesons are very similar, 

and that for heavy quarks the QQ interaction is half as strong as the 
QQ interaction. We get bbqq (1=0, J=l) bound by about 100 MeV 
and ccqq unbound. 
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1 Introduction 



The constituent quark model has been rather successful in describing the 
properties of individual hadrons |2|, [J. The extrapolation to two-hadron 
systems is, however, still rather uncertain and the predictions of various 
quark models sometimes differ dramatically. One source of discrepancy is 
our limited knowledge about the "correct" effective quark-quark interaction, 
and another source are different choices of model parameters designed to fit 
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different subsets of hadrons. Even the calculation of the nucleon-nucleon 
interaction in the model with two three-quark clusters is quite unreliable. 

Much can be learned by studying a simpler two-hadron system: two heavy 
mesons. The lowest states should be very narrow (longlived) since they can 
decay only weakly. Although their width is only ~ 1CT 3 eV (corresponding to 
a lifetime of ~1 ps) they are difficult to detect because of a low production 
cross section. Nevertheless, they are interesting theoreticaly, to confront 
different models. The detailed calculations in the literature |4], || rely on 
particular quark models, therefore we attempt an almost model-independent 
phenomenological estimate. 

In Sect. 2 we present our phenomenological estimate of the binding energy 
in order to support the expectation that some heavy dimesons may be bound. 
The estimate is mainly based on the assumption that the wave functions 
of the two light quarks around the heavy quark in A c , A& and around the 
antidiquark in the ccqq and bbqq dimesons are very similar. This assumpton 
implies that the heavy antidiquark in a colour triplet state acts just like a 
very heavy quark and that the 1/m corrections are neglected. 

This assumption is supported by the experimental observation that the 
masses of many meson or baryon pairs with c quark replaced by b quark differ 
by the same amount (essentially by the quark mass difference): B — D = 
3341 MeV, B a - D s = 3328 MeV, A b - A c = 3340 MeV. Q 

The second assumption is that for heavy quarks the QQ interaction is 
half as strong as the QQ interaction. This is true for the one-gluon-exchange 
interacton and popular for the confining interaction. In potential models with 
2-body interaction it is difficult to avoid the A ■ A factor in the interaction 
without closing all open channels for white clusters. One should note also 
that for heavy quarks the one-pion-exchange potential does not apply. 

In Sect. 3 we show that several refinements may give corrections of ±20 
MeV but they do not change the qualitative features and conclusions. 

We comment the results in the Conclusion where we explain why some 
previous model calculations are in agreement with our phenomenological re- 
sults and why some other calculations which disagree are irrelevant. 



1 In this paper we shall denote the masses of particles just by their names, and the tilde 
denotes a hyperfine average (e.g. D = jD + jD*). 
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2 The phenomenological relation for the bind- 
ing energy of dimesons 

The results in this section, Eqs. (|5]) and ([8)), are the main message of our pa- 
per. In next section we only scrutinize the assumptions and neglected terms 
and estimate the minor corrections. While in the present section the calcula- 
tions are very simple, essentially interpolations between experimental hadron 
masses, the corrections in next section rely strongly on detailed constituent 
quark model calculations 0. 

We derive phenomenological relations between dimesons, mesons and 
baryons by making two assumptions, (i) the equality of the wavefunction 
of two light quarks in the dimeson and in a heavy baryon, and (ii) the sim- 
ilarity between the subsystem of two heavy quarks in the dimeson and the 
quark- ant iquark system in a heavy meson. 

The first assumption implies that the heavy antidiquark in a colour triplet 
state acts just like a very heavy quark and that the 1/m corrections are ne- 
glected; however, there is no prejudice for the light-light subsystem in favour 
of the one-gluon-exchange or one-Goldstone-boson-exchange interaction. 

The second assumption implies that effective interaction between two 
heavy quarks in colour triplet state is half as strong as the effective quark- 
antiquark interaction in colour singlet state; therefore the binding energies 
of diquarks are weaker than those of mesons, and the distances increase. 
We do not prejudice the radial dependence of the effective interaction (the 
flavour-dependent spin-spin interaction is treated later as a perturbation). 

We call the u and d quarks q and the dimesons (bbqq) = T bb , (ccqq) = T cc 
(the symbol T reminds of the terminology tetraquarks) . 

The binding energy E bb of a quark and antiquark in a meson is a function 
of the reduced mass only, e.g. T = 6 + b + E bb , E bb = F{m = b/2). For the 
diquark bb the Schrodinger equation is similar as for the 66 meson with twice 
weaker interaction 



p2 + V bb 



2(6/2) 



" + v bb 



2(6/4) 



if) = E bh iP (1) 



To get the similarity, we have mimicked the kinetic energy with half smaller 
reduced mass. The binding energy is then E bb = \F{b/A). 
Now we compare the following hadrons 



T bb = 2b + 2q + E bb + E, 
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T = 2b + E bb 
A b = b + 2q + E qqQ , 



(2) 



where E qq q m E qq(bb ^ m E qqb is the potential plus kinetic energy contribution 
of the two light quarks in the field of a "heavy quark". The quantity E qq q 
is taken from the experimental data on A c and A b and we do not need to 
calculate it. 

We obtain the phenomenological relation 

T bb = A b + IT + 6E bbl 5E bb = ±[F(6/4) - F(b/2)}. (3) 

An analogous comparison gives for the charm-charm dimeson 

T cc = A c + + 5E CC , 5E CC = I[F(c/4) - F(c/2)}. (4) 

The binding of the (J = 0, J = 1) dimesons is expressed with respect 
to the corresponding thresholds. In order to be unambiguous we give in 
Appendix all experimental data and all references to model parameters which 
we need now or later. 

AT bb = k b + \T -B-B* + 5E bb = -250MeV + 8E bb , 

AT CC = A c + I J/tjj - D - D* + 6E CC = -42 MeV + 5E CC (5) 

Now comes an important idea how to obtain phenomenologically the "cor- 
rections" 5E. In Fig. (ED) we plot the phenomenological binding energies ob- 
tained from experimental meson masses and from a popular sets of quark 
masses, as a function of the reduced masses of the 66, cc, 6s, cs, bq mesons. We 
chose the quark masses proposed by Bhaduri et sl.(BD) and by Silvestre- 
Brac (AL1) 0], collected in Appendix; because of similarity between the two 
sets only the Fig.(|l|) for the set (BD) is presented. The results are only 
weakly dependent on quark masses. We then interpolate in order to obtain 
the values of F for the scaled reduced masses of diquarks. The idea to in- 
terpolate between experimental meson masses has been used before 0, but 
their construction of dimesons is quite different: a light diquark and a heavy 
antidiquark bound as a meson; however, the assumption of a light diquark is 
less justified than our construction with full freedom for the two light quarks. 

For the interpolation we use the following data 

\E bb = |T - 6 = \F(\b) = -529 MeV (BD), -497 MeV (ALl), 
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Figure 1: Interpolation for the bb binding energy - BD parameters (The 
curve for AL1 parameters is similar.) 

\$-c= \F{\c) = -336 MeV (BD), -302 MeV (AL1), 
\(D a -c-s) 

\F((c~ l + s" 1 )" 1 ) = -197MeV (BD), -168 MeV (AL1)(6) 

Here we used the hyperfine averages ip = ~J/i/) + jT] c and similarly for D s ; 
for T we take the S = 1 state since t] b is unknown and the splitting is anyway 
small. 

For the interpolation between T and ip we plot in Fig. (|TJ) curves of the 
form crm -1 / 3 + r. Such interpolation corresponds to the scaling for the linear 
quark- ant iquark potential; it also gives a rather straight line for phenomeno- 
logical binding energies if plotted against the abscissa m~ 1//3 . We get 

BD: \F(\b) = 1^(1315 MeV) = -407 MeV. 
All: \F(\b) = |F(1307MeV) = -375 MeV. (7) 

There is some uncertainty in the interpolation. The extreme choices of linear 
interpolation of F versus 1/m (or m) give 15 MeV stronger (27 MeV weaker) 
binding, respectively, suggesting an error ±15 MeV. The binding energy 
AT bb is then 

BD : SE bb = +122 ±15 MeV, AT bb = -128 ± 15 MeV 




2 cc 

If - — 

2 cs ~ 
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ALl : 5E bb = +122±15MeV, AT 6b = —128 ± 15 MeV (8) 

These values are very close to the result AT bb = —131 MeV of a detaild 
4-body calculation with BD interaction [|J]. 

For T cc it is even easier to interpolate since the reduced mass c/4 happens 
to lie very close to the reduced mass of the D s = cs meson so that 

F{\c) = F(467 MeV) » F^c" 1 + s -1 ) -1 ) = F(454 MeV) (9) 

for BD and similar for ALl. However, the choice of the strange quark mass 
brings some slight additional model-dependence. 
We get 

BD: 5E CC = +139 MeV, AT CC = +97 MeV 
ALl : 5E bb = +134 MeV, AT 66 = +92MeV (10) 

Therefore we expect rather reliably that the DD* system T cc is unbound. 

It is amusing to perform an analytic calculation of SE by mimicking the 
a/r + k + Ar potential by a logarythmic potential Uln(r/r ). The ground 
state energy of mesons is 

E = F(m) = U (e - \ ln[2mr^/^ 2 ]) (11) 

where m is the reduced mass and e = 1.0443 is obained by the numerical 
solution of the Schrodinger equation in dimensionless form. Then we obtain 
for any quark Q = c,b, ... 

5E QQ = |[F(Q/4) -F{Q/2)} = ±fU = 127 MeV (12) 

which is surprisingly close to the above phenomenological estimates for the 
bb and cc diquarks. The strength U = 733 MeV is taken from ref . M . 

3 Refinements 

In this section we make several refinements and corrections in order to test 
our assumptions and approximations. It turns out that the refinements may 
change the results only by 10-30 MeV so that our qualitative conclusions 
remain valid. 
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3.1 Sensitivity to quark masses 

We repeat the calculation with a rather different choice of masses. We explore 
smaller masses. We propose the following phenomenological choice in which 
we assume that the wave functions of both partners are the same. Then mass 
differences of mesons correspond to mass difference of quarks, and the ratio 
of hyperfine splits corresponds to the mass ratio of quarks. 

b- c ^ B - D = 3341 MeV, b = 4941 MeV 

b/c « ° ~° = 3.087, c = 1600 MeV 
1 B* - B 

The results (using the linear interpolation of E versus m~ 1//3 ) are 

E bb = -128 MeV, 5E bb = 83 MeV, AT bb = -167 MeV 

As we see, smaller quark masses yield a somewhat stronger binding. The rea- 
son is indirect, the separation of meson masses in quark masses plus binding 
energy is to some extent arbitrary, smaller quark masses mean less negative 
binding energy and therefore a smaller SE. 

For the strange quark such suggestions are very uncertain since the as- 
sumption of equal wave functions is not justified. Nevertheless, we take the 
suggestions as a guidance (we take for the light quark q = 320 MeV.) 

c- s ~ D - K = 1179 MeV, s ~ 421 MeV 

s/ c ~ — — — = 0.355, s ~ 568 MeV 
' K* - K 

s - q ~ D s - D = 103 MeV, s ~ 420MeV 
s - q ~ B s - B = 90 MeV, s ~ AlOMeV 

and choose for the phenomenological calculation s = 420 MeV. Since the 
reduced mass for the cc diquark m = c/4 = 400 MeV is close to the reduced 
masses for cs (333 MeV) and bs (387 MeV) we interpolate linearly and we 
obtain 

E cc = +19 MeV, 5E CC = 85 MeV, AT CC = +43 MeV 

The DD* dimeson is unbound also for these smaller masses. We conclude 
that there is some model-dependence on masses, but not too strong. 
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3.2 The spin-spin interaction 



Since the bb diquark is colour antisymmetric and orbital symetric it must 
be spin symmetric (S=l) and is analogous to T. On the other hand, to get 
F{b/A) we had to interpolate between F{b/2) and F(c/2). Since the spin-spin 
interaction is flavour dependent it scales differently with masses as the spin- 
independent part of interaction and we have interpolated between hyperfine 
averages. Therefore Fs=i(c/2) corresponding to the bb system should be 
lifted by \ ■ \{J/ijj - 7] c ) (c/b) = 5 MeV (for BD) and F(6/4) even less, which 
is negligible. (Since rjj, is still unknown, we have actually used T instead of 
T and the lifting is still smaller). 

For the cc diquark the shift is sligtly larger. In the hyperfine average, 
E cc is close to \E cS . The splitting D* — D s = 144 MeV suggests the splitting 
between cc states \{D* S — D s )(s/c) = 23 MeV. The 5=1 diquark will then 
be 6 MeV higher than the hyperfine average and then also the estimate for 
the T cc energy comes 6 MeV higher. But T cc is anyway unbound by almost 
100 MeV. 

3.3 The centre-of-mass motion 

In our phenomenological estimation we took the kinetic energy of the relative 
motion of the light diquark against the heavy diquark in T bb to be equal to 
the kinetic energy of the light diquark against the b quark in A&. To evaluate 
the correction to the kinetic energy, we repeated the calculations carried out 
in || where the BD parameter set was used. 

For the orbital part of the Tbb wave function we took a simple ansatz 
where the wave functions of heavy diquark, light diquark and their relative 
motion are Gaussians with widths bi which are chosen so that the the energy 
of is minimal 

R = ex P (-(r 1 2 /26? + r 2 2 /2foi + r 3 2 /26?)), 

b x = 0.59 fm, b 2 = 0.23 fm, b 3 = 0.60 fm. 

Here r±, r 2 and r 3 are the distances between the heavy quarks, between 
the light quarks and between the centres-of-masses of the two clusters, re- 
spectively. This wave function is symmetric in permutation of the two light 
quarks and in permutation of the two heavy quarks so that the spin-colour 
part of the function for T&& in the state with total isospin, spin and parity 
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ISP = 10 + must have the form: 



S qq C qq S bb C bb — (0) qq (3) gq (l) bb (3) bb 

Since we expect that this is the most important configuration we neglect 
the contribution of other configurations in the calculation of the correction 
to the kinetic energy. We also expect that a more elaborate wave function 
of the type which gives ~ 30MeV smaller mass of T bb would not change 
this correction by more than a few MeV. 

The calculated kinetic energy of relative motion of (bb) (qq) is 253 

MeV. This kinetic energy should be smaller than the corresponding kinetic 
energy of b (qq) in A b used in our estimation by a factor 

1/20 + 1/6 

' 1.06. 



l/2g + l/26 

Therefore the mass of T bb should be for 15MeV smaller than the phenomeno- 
logicaly estimated mass. 

We also confirmed that the width of the spatial wave function of the 
heavy diquark in T bb (calculated in the approximation with one Gaussian) 
is the same as in an isolated system of two b quarks. This supports our 
assumption that the presence of light quarks in dimeson does not change the 
binding energy of the heavy diquark (the change is less than 1 MeV). 



3.4 The finite size of the heavy diquark 

We again assume that the most important configuration in T bb is one with 
antitriplet colour function for the light diquark and triplet colour function 
for heavy diquark. This makes T bb very similar to the baryon A b which we 
use in our phenomenological estimation. Here we calculate the correction 
to the heavy-light interaction energy because of the finite size of the heavy 
diquark. If we make the limit 6 2 — > in the wave function defined in the 
previous subsection the binding energy is by 18 MeV less negative (for BD 
parameters, and a similar value for AL1). So this correction has the opposite 
sign as the correction due to the smaller kinetic energy and they almost 
cancel. 
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3.5 The colour configuration mixing 

In Tbt, the light diquark and the heavy diquark can also be in colour (anti) sextet 
state. Such additional configurations are not present in A& and can lower the 
energy of T^. This effect was calculated in Born-Oppenheimer approxima- 
tion since such approximation is very suitable for studing simultaneously the 
significance of sextet configurations as well as the tendency for two-cluster 
configurations, (next subsection). 

The wave function of is taken as a linear combination of seven config- 
urations $j listed in table [I] (for we choose an orthogonalized basis) 

*(r) = 5>,(r) $.,(!•). 

i 

We solved accurately seven coupled differential equations 

" ^ Q(r) + £ Vy(r)c,-(r) = E Ci (v) (13) 

where m^B is the reduced mass of the two B = bq clusters and the effective 
Born-Oppenheimer potential Vy(r) contains the kinetic energies of the light 
quarks as well as all interaction energies. The calculation is done in the spirit 
of the generator coordinate method (or resonating group method) except that 
Vij(r) is calculated with fixed heavy quarks. 

With the BD parameter set the binding energy is -70MeV if all seven 
configuration from table [I] are used, and -45MeV if only the first three con- 
figurations (with colour triplet wave functions) are considered. We find it 
plausible that colour sextet configurations, which were neglected in our phe- 
nomenological estimation reduce the mass of T&& by ~ 25 MeV also in an 
accurate calculation. The effect is somewhat larger in the case of T cc due to 
the smaller mass in the denominator of the spin-spin interaction which mixes 
the configurations, but the energy is still far above the threshold. 

3.6 Search for a two-cluster configuration ("molecule" 

BB*) 

At short distance, the colour triplet configurations $ 3+ ,$i + and <3>3_ give a 
Coulomb-like attraction between the two B = bq clusters while the colour 
sextet configurations give repulsion. At intermediate distances, on the other 
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hand, one can gain energy with a strong mixing between triplet and sextet 
configurations. It turns out, however, that in the Born-Oppenheimer wave 
function of the previous subsection, the amplitudes of colour sextet configu- 
rations are more than ten times smaller than the amplitudes of colour triplet 
configurations. This has the consequence that the maximum probability of 
the relative motion occurs at the origin. Detailed calculations || gave in 
fact no bound states of with a two-cluster ("molecular" or "covalent") 
structure. 

In Fig. H] we show the "adiabatic' Born-Oppenheimer potential, where 
at each distance r separately the matrix has been diagonalized. This 
is not an accurate approximation, but it demonstrates well the point that 
no second minimum appears in the BB* effective potential. Of course, the 
actual accurate calculations of eigenenergies were performed with the full set 
of equations flI3p. 




Figure 2: The adiabatic Born-Oppenheimer potential for the BB* system 



4 Conclusion 

Assuming the equality of the wave function of two light quarks in the dimeson 
and in A;, as well as the similarity between the subsystem of two heavy quarks 
in the BB* dimeson and the quark-antiquark system in botomium (with 
scaled interaction strength) we predict a binding of the = BB* dimeson 
of about -100 MeV. Analogous comparison with J/tp and A c predicts the 
dimeson T cc = DD* to be unbound. 

The binding energy of the BB* dimeson of about -100 MeV and no bind- 
ing for DD* agrees well with some previous calculations |5j in the con- 
situent quark model, since their model parameters reproduce well T, J/ip, 
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A b and A c . On the other hand, model calculations which used unrealisti- 
cally strong pion-exchange interacton yielded very different predictions for 
the dimesons (many bound states); they were of course irrelevant since they 
would strongly overbind A& and A c too. 

Experimental verification of our predictions would be a challenge since it 
would support the u Vqq = \Vqq rule" 
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A Input data and model parameters 

In order to be unambiguous we list in table |2| the baryon and meson masses 
which we use as input for our phenomenological estimate. 

We also list in tables § and |] the model parameters of Bhaduri et al. 
IJ and of Silvestre-Brac and Semay |2| which we need in Sect. 3 for our 
refinements. To unify the notation we write down the form of the effective 
quark-quark interaction 
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\ Xj j a s a s 2n{hcf 3 \ 
• — { h k + Ar H — <7j • Gj o (r) } 

exp(— r/r ) exp(— r 2 /r(, 2 ) , / 2mirrij 
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Table 1: Seven configurations for T bb with ISP = 10 + . The spatial wave 
function of T bb is constructed from two delta functions for the two heavy 
quarks and two Gaussians for the two light quarks. The lower index ± 
characterises the per mutational symmetry of two identical particles. 



M(b) \ , ( , /0 , m(b) 1,1,3 / (r,±r/2) 2N 

1 -d(r b ±r/2), ,)/ > = (—=-) 2 exp 1 



N(b) I n (b) j y ^k> F V 2k2 



0i+ = 75 [m(qi)m(q 2 ) + n(q 1 )n(q 2 )\ 

02+ = ^[m(q 1 )m(q 2 ) -n{q 1 )n{q 2 )\ 

03+ = 775 [m (qi) n (q 2 ) ± n (qj m (q 2 )\ 

V>± = j= 2 [M(b l )N(b 2 )±N(b l )M(b 2 )} 
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Table 2: Masses of some baryons and mesons in MeV. Third and sixth 
columns give the hyperfine average of the vector meson and the pseudoscalar 
meson above it 





5624 




A c 


2284.9 










Vc 


2979.8 




T 


9460.4 




JH> 


3096.9 


3067.6 


B s 


5369 




D s 


1969 




B* s 


5416 


5404 


D* 


2112 


2076 


B 


5279 




D 


1867 




B* 


5325 


5314 


D* 


2008 


1973 


K 


495.7 










K* 


893.9 


794.3 









Table 3: Choice of quark masses 



Author 




Bhaduri 


Silvestre- 
-Brac 


m u4 = q 


MeV 


337 


315 


m s = s 


MeV 


600 


577 


m c = c 


MeV 


1870 


1836 


772b = b 


MeV 


5259 


5227 



Table 4: Choice of parameters for the quark-quark interaction 



Author 




Bhaduri 
(BD) 


Silvestre- 
-Brac (ALl) 


a s 


MeV fm 


77.0 


75.0 


K 


MeV 


-685.1 


-624.1 


A 


MeV/fm 


705.7 


628.4 


a s 


MeV fm 


462.0 


257,37 


ro 


fm 


0.4545 




A 


GeV- 1 




1.6553 


V 






0.2204 
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